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I. INTRODUCTION 



II. THREE GLUONS 



The amplitudes for the production of charmonia states 
in hadron-hadron and photon-hadron collisions are usu- 
ally calculated within the framework of non-relativistic 
quantum chromodynamics (NRQCD). Several pertur- 
bative QCD reactions are required among them being 
g + g — > g + charmonia and 7 + 5 — > g + charmonia, where 
g represents a gluon. The latter can either be color singlet 
or color octet states. Specific results have been presented 
in [l|, [U, [1], [1] and [f|, among others. However a close 
examination of these papers reveals inconsistencies be- 
tween the published results. Also while the individual 
helicity amplitudes are available in the color singlet case 
we could not find the corresponding results for the color 
octet case. Therefore we have calculated the amplitudes 
by helicity methods and present our results below. For 
the benefit of the reader we also give some details of the 
calculation. 

We used the helicty method described in the book by 
Gastmans and Wu [|| (see also [l|) to calculate processes 
where three gluons or two gluons and a photon form char- 
monium. Like Gastmans and Wu we projected out the 
lowest angular momenta states of the heavy quark pair, 
namely 1 Sq, 3 Si, 1 P\, 3 Pq, 3 Pi and 3 P2, using appropri- 
ate projection operators (see pj). We then flipped one 
of the gluons from incoming to outgoing and with these 
squared matrix elements calculated the polarized and un- 
polarized differential cross sections. 
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Gastmans and Wu have presented results for the dif- 
ferential cross section for the production of a color singlet 
heavy quark pair in angular momentum states 2S+1 Lj. 
They begin with the reaction with three incoming glu- 
ons where the momenta and colors of the particles are 
labelled as 

g(k U a)+g{lS2, b)+g(k 3 , c) -> q{p/2 + q) + q{p/2-q) . (1) 

There are six Feynman diagrams where the three gluons 
couple directly to the heavy quark line and six diagrams 
where two gluons couple to the heavy quark line. 

There are eight helicity matrix elements which are la- 
belled by assigning either a + or a - to each gluon and 
which are related by CP conjugation and crossing. All 
eight can be derived from two, called |M(+, +, +)| 2 and 
\M(+, +, -)| 2 . We will list them below. 

The gluon helicities for the 2S+1 Lj (+,+,+) combina- 
tion are 

£ = N\M 2 U^~ 75) +#yyi(l+7 5 )] (2a) 

# = iV[^i(l-75) + ^ 2 (l + 75)] (2b) 

# = A[^ 2 (l- 75 ) + ^^ 3 (l + 75 )], (2c) 

while those for the +1 Lj (+, +, — ) combination are 

£ = N[^J 3 (l--y 5 ) + ^J 1 (l + y 5 )] (3a) 
4 = -Ar[^ 3 (l- 75 ) + ^^ 2 (l + 7g )] ( 3 b) 
& = iV[^^ 2 (l+ 75 ) + ^^ 3 (l- 75 )], (3c) 

where N = [(ki ■ k 2 )(k 2 ■ k 3 )(k 3 ■ fc x )]-*/4 is a normal- 
ization factor. In principle there should be extra terms 
in these expressions but they do not change the answers 
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in this reaction due to the symmetric choice of variables. 
The helicity amplitudes are functions of the invariants 
s = (ki+k 2 ) 2 , t = (k 2 + k^) 2 , u = (fc 3 + /ci) 2 and the mass 
of the pair M ~ 2 m where m is the heavy quark mass. 
Note that s + t + u = M 2 , N = (2stu)^^ and the color 
singlet projection operator is %/\/3- Also they depend 
on two parameters Rq and R[ which are the S-state wave 
function and the derivative of the P-state wave function 
evaluated at the origin. The former is defined in terms 
of the leptonic decay width 



R 2 = M 2 r( 3 5i -» e + e-)/4a 2 Q 2 f , 



(4) 



with a w 1/137 the fine structure constant and Qf is the 
fractional charge of the quarks. R[ is determined from a 
fit to the charmonium potential and has the value 



R[ /M 2 w 0.006 (GEV) E 



(5) 



We actually need the differential cross section for the 
reaction 



Sij I \/3 in the color singlet case is replaced by V2T^ . Also 
the color octet amplitudes cannot be determined from 
decay processes so they are fit to quarkonium produc- 
tion differential cross sections in proton-proton, proton- 
antiproton and photo-hadron collisions. 

We compare our results with those in previous papers. 
The differential cross section for unpolarized reactions 
such as P+P — > cc+X contains the sum of the squares of 
the helicity amplitudes, |M(+,+;+)| 2 + |A/(+,+:-)| 2 
+ |M(+,-;-)| 2 + |M(-,+;-)| 2 . 

The color singlet case results are given by @ and [lj, 
which we refer to as GW and GTW respectively. The 
color octet results are available in the Appendix of Cho 
and Leibovich which we refer to as CL. The differen- 
tial cross sections for longitudinally polarized collisions 
contain the differences \M(+, + ; +)| 2 + \M(+, +; -)| 2 - 
|M(+,-;-)| 2 - |Af(-,+;-)| 2 , and are listed for both 
the color singlet and the color octet cases in the paper 
of Klasen, Kniehl, Mihaila and Steinhauser H, which we 
refer to as KKMS. 



/(fci, a)+g{k 2 , b) -> q(p/2+q)+q(p/2-q)+g(k 3 , c) , (6) 



where the invariants are now s = (ki+k 2 ) 2 , t = (k 2 — k^) 2 , 
u = (ki - k 3 ) 2 . 

The squares of the matrix elements for the reaction ^ 
follow from those in reaction fl} by crossing k% — > — k$ 
and flipping the helicity of the third gluon. 

They are denoted by M|(+,+;+)| 2 , |M(+,+;-)| 2 , 
|M(+,-;-)| 2 and |M(-,+;-)| 2 . Note that these are 
equal to M|(-,-;-)| 2 , |M(-,-;+)| 2 , |Af(- +;+)| 2 
and |M(+, — ; +)| 2 respectively by CP conjugation. How- 
ever the kinematic variables require permutations to re- 
flect the crossing of gluon number three. These relations 
are 



\M(- 
\M(- 



;-F = 



\M(+,--)\ 2 = 



iM(+,+,-)r 



|M(- +;-)| 2 = |M(+,+,-)r 



(7a) 
(7b) 
(7c) 

(7d) 



It is very convenient to use completely symmetric vari- 
ables which are then invariant under any crossing trans- 
formations. Hence we express several results in terms of 
the variables M 2 = s + t + u, P = st + tu + us and 
Q = stu, which are invariant under s <-> t and s «-» u. 
The denominators of the helicity amplitudes are written 
in these variables while the numerators contain terms in 
s. Therefore the crossing simply involves changing s — * t 
and s — > u in the numerators of our expressions. 

We have also calculated the corresponding amplitudes 
for the production of a color octet heavy quark pair which 
requires four additional Feynman diagrams for processes 
where only one gluon couples to the heavy quark pair. 
These contain three gluon and four gluon couplings. The 
color octet projection operator is required so the factor 



A. Matrix Elements Squared 

We now list the results for the squares of the color sin- 
glet matrix elements when the heavy quark pair (with 
mass M) is in the appropriate angular momentum state. 
However for convenience we rename i? 2 , = (R^Sq] ) = 



R[ 3 S[ 1} }) and R' 2 = 
\ - ' R[ 3 P 2 {1) ] 



(i)i 



R?P ( 



Mi 



), where the final superscript in- 



dicates the color singlet. 



For 1 Sq we find 

|M(+,+,+)| 2 = 
|M(+,+,-)| 2 = 



Color Singlet 



lc(l)l 



M 8p2 



irM Q(Q - M 2 P) 2 



s 4 P 2 



ttM 



Q(Q-M 2 P) 2 



(8a) 
(8b) 



where the color states of the gluons have been summed 
over. These results agree with the squares of (8.29) and 
(8.40) in GW. 
For 3 Si we find 



\M(+,+,+)\' 

iM(+,+,-)r 



= (9a) 
_ 160g 6 ( R^S^} ) M 2 s 2 {s - M 2 ) 2 



9nM 



(Q - M 2 P) 2 



-,(9b) 



where the color states of the gluons have been summed 
over. The polarization of the spin one charmonium state 
has also been summed over. The second result agrees 
with the (8.50) in GW after correcting an obvious typo 
that the (t - M 2 ) should read (t - M 2 ) 2 . 
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For 1 Pi we find, after summing over colors and polar- 
izations, 



(10a) 



M 2 s 



2 Jl 



640. 9 6 ( P^P^]) 
3ttM 3 
M 10 (-M 2 P + 5Q) 
X (Q - M 2 P) 3 

640. 9 6 (i?[ 1 P 1 (1) ] 

3ttM 3 (Q - M 2 P) 3 
x [3M 4 Q - M 6 P + 2Qs 2 ] . (10b) 

Here we agree with the results (8.55) and (8.57) in GW. 

For 3 P we find, after summing over colors and polar- 
izations, 

.r3 p(i)i 



|M(+,+,+)| 2 = 



|M(+,+,-)| 2 = 



9M 8 P 2 (Q -M 2 Pf 



m+,+,-)\ 



Q(Q - M 2 Py 

2 = 64g 6 (P[ 3 P (1) ]) {s-M 2 ) 2 
nM 3 



(11a) 



Q(Q - M 2 P) 4 
Q 2 - s 2 Q(s - 3M 2 ) + 3PM 2 s 3 ] 2 . (lib) 



Here we agree with the results in (8.59) in GW. 

For 3 Pi we find, after summing over colors and polar- 
izations, 

|M(+,+,+)| 2 







(12a) 



192 5 6 (P[ 3 P 1 (1) ]) (s-M 2 ) 2 s 



\M(+,+,-)\ 



2\2 2 



irM 3 (Q - M 2 P) 4 

2Q(5M 4 P - M 8 + P 2 
(4P - 2sM 2 + As 2 - M 4 ){s - M 2 ) 2 ) 
Q 2 {lhM 2 - 8s) - 4M 2 P 3 + M 6 P 2 ] , (12b) 



which agrees with (8.63) in GW. 

For 3 P 2 we find, after summing over colors and polar- 
izations, 



\M(+,+,+)\ 
\M(+,+ 





2 _ 64 5 6 (P[ 3 P 2 (1) ]) 



n irM 3 Q(Q-M 2 P) 4 

12M 8 P 4 (3s — M 2 ){s — M 2 ) 

-l2M 4 P 5 s(s - 3M 2 ) + 2P 2 Q 3 (s - 11M 2 ) 
-3M 6 P 3 Q(s - M 2 )(25s - 8M 2 ) 
+12M 2 P 4 Q(s 2 - 4Af 2 s - 3M 4 ) 
+M 4 P 2 Q 2 (8s 2 + 9M 2 s - 15M 4 ) 
-2P 3 Q 2 (s 2 - 5M 2 s - 30M 4 ) 
+M 2 PQ 3 (29s 2 - 51Af 2 s + 18M 4 ) 

-M 2 g 4 (9s - 11M 2 )] , (13b) 
which agrees with (8.70) in GW. 



2. Color Octet 



Now we present the corresponding results for the color 
octet projections. These results do not seem to be avail- 
able in the literature. We have only found expressions for 
the differential cross sections which we will compare to 
ours later on. We give these results since we need the dif- 
ferences between the helicity combinations to check the 
octet longitudinally polarized differential cross sections. 
The constants from the wave functions are now simply 
renamed as P 2 — > (R ^Sq 8 ^] ) etc., since there are other 
definitions in the literature. We will present the relations 
between the definitions later on. 



For 1 S'o we find 



\M(+,+,+)\ 



2 _ 40g 6 (P ^S^} ) M 8 (P 2 - M 2 Q) 



ttM Q(Q - M 2 P) 2 

(14a) 



2 40 3 6 (P[ 1 ^ 8) ]> 



ttM Q(Q-M 2 P) 2 
x PQ + s 3 (s - M 2 ) 2 - s 2 q] . (14b) 



For 3 Si we find 



\M(+,+,+)\ 2 = 
\M(+,+,-)\ 2 = 



For 4 Pi we find 



(15a) 

16g 6 (R[ 3 s{ 8) }) s 2 (s-M 2 ) 2 

3irM M 2 (Q - M 2 P) 2 
x(19M 4 - 27P) . (15b) 



(13a) l^(+,+,+)| 2 = 



32.g 6 (P[ 1 P| 



lp(8)l 



M 6 



ttM 3 Q(Q - M 2 P) 3 
2\7M 4 Q 2 - 54PQ 2 + 43M e PQ 

-27M 2 P 2 Q - 27M 4 P 3 

, |2 32g 6 (P[ 1 P 1 (8) ]) s 3 
|M(+,+,-)| 2 = - y ! 



(16a) 



ttM 3 Q(Q - M 2 P) 3 

Qs 2 (t + w)(-174u 2 + 26tu - \7At 2 ) 

+Qs(-98u 4 - 278tu 3 - 468iV - 278t 3 u - 98t 4 ) 
+Q(t + u)(-38u 4 - 82tu 3 - 169iV - 82t 3 u - 38t 4 ) 
+s 4 {-27u 4 - I52tu 3 + lOiV - lh2t 3 u - 27t 4 ) 



+t 2 u 2 (t + u)(-38u 3 - 60tu 2 - 60t 2 u - 38t 3 ) 
+s 5 (t + u){-27u 2 -lltu-27t 2 ) . 



(16b) 
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For 3 P we find 



2 160g 6 (P[ 3 P (8) ] 



\M(+,+, 



irM 3 

9M 8 (P 2 - M 2 Q)(Q - M 2 P) 2 

Q(Q - M 2 py 

160g 6 (P[ 3 P (8) ]) (s-M 2 ) 2 



(17a) 



ttM 3 Q(Q - M 2 P) 4 

x Q 4 + 9s 8 M 4 (s-M 2 ) 2 

+Qs 5 M 2 (6s 3 - 6M 6 + 33sM 4 - 42s 2 M 2 ) 
+Q 2 s 2 {Us 2 M 4 + AM 8 - 18s 3 M 2 + s 4 ) 

+g 3 s(-2(s - M 2 ) 2 + 9sM 2 )] . (17b) 



For 3 Pi we find 
|M(+,+,+)| 2 = 



(18a) 



\M(+,+,-)\ 



2 _ geoff^ppp^]) 



ttM 3 

(s-M 2 ) 2 (P + s 2 -sM 2 ) 



Q 3 (s - 2M 2 ) 



Q(Q-M 2 Py 
+s 5 M 2 (M 8 - 4sM 6 + 7s 2 M 4 - 6s 3 M 2 + 2s 4 ) 
+Qs 3 (M 8 - 4sM 6 + lls 2 M 4 - 10s 3 M 2 + s 4 ) 



+Q 2 s(M 6 + 7s 2 M 2 - 2s 3 ) 



(18b) 



For 3 P 2 we find 

|M(+,+,+)| 2 = (19a) 

_ 2 320g6(P[ 3 Pf]) (s-M 2 ) 2 
1 1 ' ' Jl ' ttM 3 Qs 4 {Q-M 2 P) 4 

x 6s 8 M 4 (s - M 2 ) 6 + Qs 6 M 2 (18M 12 - 114sM 10 

+285s 2 M 8 - 354s 3 M 6 + 225s 4 M 4 - 66s 5 M 2 + 6s 6 ) 
+Q 2 s 4 (24M 12 - 132sM 10 + 313s 2 M 8 - 336s 3 M 6 
+ 161s 4 M 4 - 30s 5 M 2 + s 6 ) + Q 3 s 2 (18M 10 - 78sM 8 
+ 141s 2 Af 6 - 110s 3 M 4 + 25s 4 M 2 - 2s 5 ) 

+Q 4 (s 4 - 6M 2 (s - M 2 ) 3 + 6sAf 4 (s - M 2 ))} . (19b) 



B. Unpolarized Differential Cross Sections 

These follow from the sum of the squares of the he- 
licity matrix elements |M(+,+;+)| 2 + |M(+,+;-)| 2 + 
|M(+,-;-)| 2 + |M(-,+;-)| 2 with the substitutions 
s — > t and s — > u as described above. However to sum 
over all polarization states we have to multiply by 2 to 
include the CP conjugates. Then one adds the average 
over the initial gluon colors and polarizations (1/256) and 
multiplies by an overall factor of l/(167rs 2 ). 



1. Color Singlet 

These results can be compared with the results in GW 
and KKMS. The latter authors give the differential cross 
sections as functions of polarization factors £ a £b in the 
form a(s, t, u) +^ a ^b b(s, t, u). The unpolarized cross sec- 
tions are obtained by setting £ a £b = 0. We call these the 
first terms and the coefficients of £a£b, which yield the 
longitudinally polarized differential cross sections, the 
second terms. Note that, due to the differences in the 
definitions of the wave functions, our comments concern 
the polynomial dependence of a{s,t,u) and b(s,t,u) on 
the invariants. However we will also identify the prefac- 
tors. This is possible because their polarized differential 
cross sections agree with ours. 

For 1 Sq we find 



da 
"dT 



Ms 2 Q(Q - M 2 P) 2 
x (P - M 4 ) 2 + 2M 2 q] , 



(20) 



which agrees with (8.46) in GW. However it docs not 
agree with the first term in (A. 16) in KKMS, who use 
the notation where (R^S^} ) = 47r(0[ 1 sf } ]). 
For 3 S*i we find 



da 
dt 



107ra 3 ( R[ 3 S[ 1] ] ) M 2 (P 2 - M 2 Q) 



9Ms 2 



(Q - M 2 py 



(21) 



which agrees with (8.52) in GW and also agrees with 
the first terms in (A. 17) in KKMS, who use the notation 
where ( P[ 3 ^ 1} ] ) = 4tt(0[ 3 S 1 1) ]) /3. 
For 1 Pi we find 



da 
dt 



40 7 ra 3 (P[ 1 P 



Mi 



M 2 



M W P 



3M 3 s 2 [Q - M 2 P) 3 L 
+M 6 P 2 + Q(5M 8 - 7M 4 P + 2P 2 ) + 4M 2 Q 2 ] , (22) 



which agrees with (8.58) in GW. It also agrees with the 
first terms in (A. 18) in KKMS, who use the notation 

where (PfP^]) = 4^(0[ 1 P 1 (1) ]) /9. 
For 3 P we find 



da 47ra 3 (P[ 3 P 
di " 



1 



2M 8 P 2 Q 



M 3 s 2 Q(Q - M 2 P) 4 
+6M 6 P 3 Q(3P - M 4 ) - 2M 2 PQ 3 {P - M 4 ) 

+P 2 (3PM 2 - Q) 2 (P - M 4 ) 2 + 6M 4 Q 4 1 , (23) 



which agrees with (8.60) in GW. It does not agree with 
the first terms in (A. 19) in KKMS, who use the notation 
where (P[ 3 P (1) ]) = 4^(O[ 3 P (1) ]) /3. 
For 3 Pi we find 



da 
~dt 



127ra 3 (P[ 3 P| 



(!)l 



P 2 



- i5M 2 g 2 



M 3 s 2 (Q - M 2 P) 4 
+M 2 P 2 (M 4 - 4P) - 2Q(M 8 - 5M 4 P - P 2 ) , (24) 
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which agrees with (8.64) in GW. It does not agree with 
the first terms in (A. 20) in KKMS, who use the notation 

where ( P[ 3 P X (1) ] > = 4 7 r(0[ 3 P 1 (1) ]) /9. 
For 3 i-2 we find 



do- 
ck 7 



47ra 3 (P[ 3 P 2 (1) ] 



1 



APs 2 Q(Q - M 2 P) 4 
12M 4 P 4 (P - AT 4 ) 2 + M 2 PQ 3 (16M 4 - 61P) 
3M 2 P 3 Q(8M 8 - M A P + 4P 2 ) + \2M 4 Q 4 



2P 2 Q 2 (7M ei - 43Af 4 P - P z ) 



(25) 



which agrees with (8.71) in GW after correcting a typo. 
They have (8M 8 -M 4 P + P 2 ) which should read (8M 8 - 
M 4 P + 4P 2 ). The expression is given correctly in their 
published paper |l[ . Also it does not agree with the first 
terms in (A. 21) in KKMS, who use the notation where 
(P[ 3 P 2 (1) ])=4^(0[ 3 P 2 (1) ])/15. 

In view of these differences we contacted the authors 
of the KKMS paper. They calculated their results with 
projection operators for the sums over the gluon polariza- 
tion states, which required the calculation of additional 
ghost diagrams. However they inadvertently presented 
the formulae (A.16), (A.19), (A.20) and (A.21) without 
the contributions from these ghost terms. They claim 
that the correct formulae are included in their fortran 
programs and that their numerical results are therefore 
correct. 



2. Color Octet 

These can be compared with the results for the squares 
of the matrix elements in the appendix of CL and with 
the first parts of the expressions in Appendix A of KKMS. 

First we find for 1 S'o 



da 
dt 



5n^{RlS^]) P 2 -M 2 Q 
Q{Q - M 2 P) 2 

- 2M 2 Q] , 



2Ms 2 



4\2 



(P - M ) 



(26) 



which agrees with (A5a) in CL. It does not agree with 
the first part of (A. 22) in KKMS, who use the notation 
where ( R^S^} ) = tt( OfS^] ) /2. 
Next, we find for 3 Si 



da na 3 (R[ 3 s[ 8) }) (P 2 
dt 



M 2 Q)(19M 4 - 27 P) 



3Ms 2 



M 2 (Q -M 2 P) 2 



,(27) 



which agrees with the sum of (A5b) plus (A5c) in CL. It 
does not agree with the first terms in (A. 23) in KKMS, 
who use the notation where ( R[ 3 s{ 8) ] ) = tt( 0[ 3 s{ 8) ] ) /6. 



The expression for 1 Pi , 



da 



27ra 3 (R[ 1 Pl 



(8)i 



1 



M 3 t 



4/n3 



179APQ 



Q(Q - M 2 P)3 
+217M 10 Q 2 - 27A/ 2 P 5 + 54M 6 P 4 - 27M 10 P 3 
+ 135PQ 3 + 103M 2 P 2 Q 2 - 212M 6 PQ 2 



-124M 8 P 2 Q + 43M 12 PQ + 27P 4 Q 



(28) 



is not given in CL. It does not agree with the first 
terms in (A. 24) in KKMS, who use the notation where 

(PpP^]) =^((9[ 1 P 1 (8) ])/18. 
Now we turn to the expression for 3 P 



dcr 
d 7 



107ra 3 (P [ 3 P^ 



(8)l 



1 



M 3 s 2 Q(Q - M 2 P) 4 
x 9M 4 P 4 {P - M 4 ) 2 + 3M W P 3 Q - 6M 2 P 5 Q 

+27M e P 4 Q + 18M 12 PQ 2 - 32M S P 2 Q 2 
-9M 14 P 2 Q - 4M 4 P 3 Q 2 + 5M 4 Q 4 + P 4 Q 2 

+ 11M 6 PQ 3 



M 2 P 2 Q 3 - 13M 10 Q 3 



(29) 



which agrees with (A5d) in CL. It does not agree with 
the first terms in (A. 25) in KKMS, who use the notation 

where ( P[ 3 P (8) ] ) = tt( O[ 3 P (8) ] ) /6. 
For 3 Pi we find 



da 
dt 



607ra 3 (P [ 3 P| 



(8)l 



1 



M 3 s 2 (Q - Af 2 P) 4 

-Af 6 P 4 - 2A/ 2 P 5 - 2M 8 P 2 Q 

-3M 6 PQ 2 - 9M 2 P 2 Q 2 + 6M 4 Q 3 



P 4 Q + M W Q 2 
- 1M 4 P 3 Q 

(30) 



which agrees with the sum of (A5e) and (A5f) in 
CL. It does not agree with the first terms in (A. 26) 
in KKMS, who use the notation where (P[ 3 Pi (8) ]) = 

^(0[ 3 Pi (8) ])/6. 
For 3 P 2 we find 



dcr 
dt 



207ra 3 (P[ 3 P 2 (8) ]) 



1 



6AT 12 P 4 - 12Af 8 P 5 



Q(Q - M 2 P) 4 

VI 4 P 6 + 11A/ 4 Q 4 

Q(-6Af 14 P 2 - 3Af 10 P 3 + 3Af 6 P 4 - 6M 2 P 5 ) 
Q 2 (24Af 12 P- 29Af 8 P 2 



41Af 4 P 3 + P 



Q 3 (-19A/ iu + 14APP - 31A/ 2 P 2 ) 



(31) 



which agrees with the sum of (A5g) plus (A5h) plus 
(A5i) in CL, after correcting an obvious typo that the 
term —Ms 2 , which multiplies the second line in (A5i), 
should read — M 2 s. It does not agree with the first 
terms in (A. 27) in KKMS, who use the notation where 
(P[ 3 P 2 (8) ] ) = tt{ 0[ 3 P 2 (8) ] ) /30. The explanation for the 
difference between our results and (A. 22) - (A. 27) in 
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KKMS is again that they inadvertently neglected to in- 
clude ghost contributions to their amplitudes. However 
they claim that they did so in their computer programs 
so their numerical results are correct. 

In view of the differences in the above results and be- 
fore contacting KKMS we recalculated the differential 
cross sections by summing over the physical polarizations 
of the external gluons using the covariant expression 



J2 e»(k,a)e»(k,a) = P^(n,k). 



(32) 



«=+,- 



ith 



P^(n, k) = -g^ + {n^K + k^n v )/n ■ k , (33) 



where satisfies n, 



pi iv = p^ n = o and n 2 = 0. One 



uses this sum for each external gluon and the answer 
for the square of the matrix elements should be indepen- 
dent of n^. This method does not require any ghosts 
and yielded the same answers we obtained above for the 
differential cross sections. 



For 1 Pi we find 



dt 



dAa _ AOna^R^P^]) 



M 2 



3M 3 s 2 (Q - M 2 P) 3 



Atu 6 



6t z u 



2„,2 



At 3 U ■ 



2t 4 



stu(2u 4 - 

+s 2 (5tu 4 + 7t 2 u 3 + 7t 3 u 2 + 5t 4 u + t 5 + u 5 ) 
+s 3 (10tw 3 + lOiV + lQt 3 u + At 4 + Au 4 ) 
+s 4 (10tw 2 + lQt 2 u + 6t 3 + 6u 3 ) 
+s 5 {Atu + At 2 + Au 2 ) + s 6 (t + u) 

(36) 



+t V + 3tV + 3t 4 u 3 + t b u 2 



This is in agreement with the second terms in (A. 18) 
in KKMS, when we make the replacement (R^P^] ) — 
4tt<0[ 1 P 1 (1) ])/9. 

For 3 P we find 



C. Polarized Differential Cross Sections 

Now we calculate the expressions |M(+,+;+)| 2 + 
|M(+,+;-)| 2 - |M(+,-;-)| 2 - |M(-,+;-)| 2 , which 
yield the longitudinally polarized differential cross sec- 
tions. 



1. Color Singlet 

We begin with the color singlet expressions. These are 
available in KKMS as the second terms, i.e., b(s,t,u), 
those terms proportional to £ a £b- The prefactors are 
identified as in the unpolarizcd differential cross sections 
given previously. We repeat them here for convenience. 

For 1 Sq we find 



dAa na 3 (R[ 4 S^}) 



s 6 - 2Q 2 



dt ' s 2 M Qs 2 (Q - M 2 P) 2 

+AQs{s - M 2 ) 2 + +s 2 M 8 - s 2 (s - M 2 ) 4 ] . (34) 

This is in agreement with the second terms in (A. 16) 
in KKMS, when we make the replacement (R^Sq 1 ^] ) = 

For 3 Si we find 



dAa 
~dT 



10to 3 (P [ 3 s{ 1] ]) M 2 Q(s 2 - P) 
9Ms 2 s{Q - M 2 P) 2 



(35) 



This is in agreement with the second terms in (A. 17) 
in KKMS, when we make the replacement (iipS^] ) = 
4tt<0[ 3 ^ 1) ])/3. 



dAcr _ 47ra 3 (i?[ 3 P (1) ] 



s 2 M 2 



dt M 3 s 2 Qs 6 (Q - M 2 Py 

x 18s 9 M 4 (s - M 2 ) 6 + 9s 10 M 6 (s - M 2 ) 4 

+Qs 7 M 2 (66M 12 - 327 sM w + 684s 2 M 8 
-762s 3 M 6 + 462s 4 M 4 - 135s 5 M 2 + 12s 6 ) 
+Q 3 s 3 (66M 10 - 260sM 8 + 374s 2 M 6 - 5s 5 



-237s 3 M 4 + 62s 4 M 2 ) + Q 5 (6sM 2 



9M 4 



+Q 4 s(18M 8 - 75sM 6 + 80s 2 M 4 - 31s 3 M 2 + 4s 4 ) 
+Q 2 s 5 (96M 12 - 422sM 10 + 750s 2 M 8 

-663s 3 M 6 + 286s 4 M 4 - 49s 5 M 2 + 2s 6 )] . (37) 



This is in agreement with the second terms in (A. 19) 
in KKMS, when we make the replacement (i?[ 3 Pg 1 ' ) ] ) = 
4^(O[ 3 P (1) ])/3. 

For 3 Pi we find 



2s 6 (-Af 8 



dAa 12na 3 (R[ 3 P^}} Q(Q + s 2 M 2 ) 
dt ~ M 3 s 2 s 5 (Q-M 2 P) A 

+5sM 6 - 9s 2 M 4 + 7s 3 M 2 - 2s 4 ) + Qs 3 (M 8 - AsM 6 
+ lls 2 M 4 - 18s 3 M 2 + 10s 4 ) + g 2 s(M 6 - 6sM 4 

+ lls 2 M 2 - 8s 3 ) + 2Q 3 (s - 2M 2 )1 . (38) 



This is in agreement with the second terms in (A. 20) 
in KKMS, when we make the replacement (PpPi 1 ^] ) = 
Att{0[ 3 p[ 1] ]) /9. 
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For 3 P 2 we find 



dAa 4to 3 (P[ 3 P 2 (1) ]) Q + s 2 M 2 
dt ~ M 3 s 2 Qs 6 (Q - M 2 P) 4 

x - 24s 9 M 4 (s - M 2 f - 12s 10 M 6 (s - Af 2 ) 4 

+Qs 7 M 4 (-48M 10 + 276sM 8 - 648s 2 Af 6 
+768s 3 Af 4 - 456s 4 A/ 2 + 108s 5 ) 
+Q 4 s(24M 8 - 63sM 6 + 34s 2 Af 4 - 5s 3 Af 2 + 8s 4 ) 
+Q 3 s 3 (15s 3 M 4 - 79sM 8 + 28s 2 M 6 - 2s 4 M 2 
-10s 5 + 48M 10 ) + Q 5 (-12M 4 + 12sM 2 - 2s 2 ) 
+Q 2 s 6 (-330sM 8 + 306s 2 M 6 - 82s 3 M 4 

-14s 4 M 2 + 4s 5 + 116M 10 )] . (39) 



This is in agreement with the second terms in (A.21) in 
KKMS, when we make the replacement (P[ 3 P 2 ]) = 
4tt(0[ 3 P 2 (1) ]) /15. Our polarized differential cross sec- 
tions agree with those in KKMS because their method 
of calculation docs not require ghost contributions. 



2. Color Octet 

These are only available in the Appendix of KKMS. 
We begin with 1 So- 

dAa _ 57raf(g[ 1 5f ) ]) 1 



dt 2s 2 M s 4 Q(Q - M 2 P) 2 

x 2s 7 M 2 (M 2 -s) 4 + s 8 Af 4 (A/ 2 -s) 2 

+Qs 5 (4M 8 - 15sAf 6 + 20s 2 M 4 - 12s 3 M 2 + 2s 4 ) 
+Q 2 s 3 (4M 6 - 12sM 4 + 14s 2 Af 2 - 5s 3 ) 

+Q 3 s(2M 4 - 5sM 2 + 4s 2 ) - Q 4 j . (40) 

This is in agreements with the second terms in (A. 22) 
in KKMS, if we make the replacement (R ^Sq 8 ^] ) = 
n{0 ?sW])/2. 
For 3 Si we find 



dAa _ 7 ra 3 (i?[ 3 5f ) ])Q(19M 4 -27P)(s 2 -P) 



dt 



3Ms 2 



M 2 s{Q - M 2 Pf 



(41) 



which agrees with the second terms in (A. 23) in KKMS, 
if we make the replacement (R [ 3 s[ 8) ] ) = tt(0 [ 3 sf } ] ) /6. 



For 1 P\ we find 

dAa 2ira 3 s (R fP™] 
dt ~~ " 



M 2 



M 3 s 2 
(27(M 2 - 



Qs 6 (Q-M 2 P) 4 



) 3 M b s li (2M 4 - 3sM 2 + 2s 2 ) 



+Qs 9 (M 2 - s)Af 4 (-621sM 6 + 864s 2 Af 4 
-567s 3 M 2 + 108s 4 + 173M 8 ) 
+Q 2 s 7 M 2 (-1395sM 8 + 2307s 2 M 6 
-1988s 3 M 4 + 621s 4 M 2 - 54s 5 + 249 M 10 ) 
+Q 3 s 5 (-1488sM 8 + 2314s 2 M 6 - 1492s 3 M 4 
+189s 4 M 2 + 249M 10 ) + Q 4 s 3 (-779sM 6 
+1379s 2 M 4 - 449s 3 M 2 + 173M 8 ) 
+Q 5 s(-162sM 4 + 373s 2 M 2 + 27s 3 + 54M 6 ) 



-Q 6 27(M 2 -s)) 



(42) 



This is in agreement with the second terms in (A. 24) 

in KKMS, if we make the replacement (R^P^}) = 

ir(0 [ipW] )/18. 
For 3 Po we find 

107ra 3 (P[ 3 P (8) ]) 



dAa 
~dT 



1 

M 3 s 2 Qs e (Q - M 2 P) 4 

9s n M 6 (2(M 2 - s) 6 + sM 2 (M 2 - s) 4 ) 

+3Qs 9 M 4 (23M 12 - 126SAP 10 + 285s 2 M 8 
-342s 3 A/ 6 + 228s 4 M 4 - 78s 5 M 2 + 10s 6 ) 
+Q 2 s 7 M 2 (117M 12 - 612sM 10 + 1285s 2 M 8 
-1358s 3 A/ 6 + 738s 4 M 4 - 184s 5 M 2 + 14s 6 ) 
+Q 3 s 5 (117M 12 - 553sM 10 + 1021s 2 A/ 8 + 2s 6 
-881s 3 M 6 + 352s 4 M 4 - 54s 5 M 2 ) 
+Q 5 s(18M 8 - 81sM 6 + 88s 2 M 4 - 33s 3 M 2 + 4s 4 ) 
+Q 4 s 3 (69M 10 - 292sM 8 + 439s 2 M 6 - 275s 3 M 4 

+68s 4 M 2 - 5s 5 ) - Q 6 (s - 3AP 2 ) 2 ] . (43) 

This is in agreement with the second terms in (A. 25) 
in KKMS, if we make the replacement (P[ 3 P (8) ]) = 

^<O[ 3 P (8) ])/6. 
For 3 Pi we find 

dAa _ 60na 3 s {R[ 3 Pi 8) }) Q 



dt 



M 3 s 2 s 5 (Q - M 2 P) 4 



s 7 M 2 (-M 10 + 5sM s - 10s 2 M 6 + lls 3 Af 4 + 2s 5 

-7s 4 M 2 ) + Qs 5 (-2M 10 + 8sA/ 8 - 14s 2 M 6 + 2s 5 
+17s 3 M 4 - 12s 4 M 2 ) + Q 4 (-s + 2Af 2 ) 
+Q 2 s 3 (-2Af 8 + 4sA/ 6 - 8s 2 M 4 + 10s 3 M 2 

+g 3 s(-A/ 6 + 3sAf 4 - 5s 2 M 2 + 4s 3 ) . 



5s 4 



(44) 

This is in agreement with the second terms in (A. 26) 
in KKMS, if we make the replacement (R[ 3 P^ 8) ]) = 
n{0[ 3 P&})/6. 
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For 3 P 2 we find 
dAa _ 20to 3 (P[ 3 P 2 (8) ] 



1 



dt APs 2 Qs 6 {Q - M 2 P) 4 

- 12s n A/ 6 (Af 2 - s) 6 ~ 6s 12 M s (M 2 - s) 4 

+Qs 9 Af 4 (-27Af 12 + 177sM w - 447s 2 A/ 8 
+567s 3 Af 6 - 378s 4 Af 4 + 120s 5 A/ 2 - 12s 6 ) 
+Q 2 s 7 M 2 (-15M 12 + 138sAf 10 - 398s 2 M 8 



+481s 3 M 6 - 255s 4 M 4 
+Q 3 s 5 (15Af 12 + 5sAf u 



47s 5 M 2 4 
- 89s 2 M 8 



2s b ) 

f 115s 3 M 6 



-35s 4 M 4 - 12s 5 Af 2 + 2s 6 ) 

+g 5 s(12M 8 - 39sM 6 + 25s 2 A/ 4 - 3s 3 A/ 2 + 4s 4 ) 
+Q 4 s 3 (27Af 10 - 55sAf 8 + 37s 2 Af 6 - 5s 3 A/ 4 

+2s 4 Af 2 - 5s 5 ) + Q 6 (6sAT 2 - s 2 - 6A/ 4 )j . (45) 

This is in agreement with the second terms in (A. 27) 

in KKMS, if we make the replacement (R [ 3 P 2 (8) ] ) = 

^(O[ 3 P 2 (8) ])/30. 

Again we agree with the KKMS results because, in 
their method of calculation, the polarized differential 
cross sections do not require ghost contributions. 



III. ONE PHOTON, TWO GLUONS 

We also calculate the helicity matrix elements for the 
reaction 

j{k 1 )+g(k 2 ,b)+g(k 3 ,c) - q(p/2 + q)+q(p/2-q) . (46) 

Here we can compare our results with those m KKMS as 
well as those in Yuan, Dong, Hao and Chao [4, which we 
refer to as YDHC and in Ko, Lee and Soy [3], which we 
refer to as KLS. 

In this case there is no t <-> u symmetry. Also we have 
to change our choice for the helicities. We use 

= 7V^ 1 ^^(lT75)-^^^i( 1 ±75)±2fc 2 .fc 3 ^ l75 ] 



(47a) 



A%^ 2 (1 ± 75 ) + #2^3(1 T 75) - 2ki.*3jf 2 



(47b) 
(47c) 

A. Matrix Elements Squared 



The squares of the matrix elements now no longer have 
all the symmetries as in the previous case, so we need the 
four helicity amplitudes | Af (+, +, +)| 2 , \M(+, +, -)| 2 , 
|Af(+,-,+)| 2 and |Af(-,+,+)| 2 . 

When we calculate the differential cross section for 

7 (fci) +g(k 2 , b) -> q(p/2 + q) +q(p/2- q) +g(fc 3 , b) . (48) 



we have to cross gluon number three. For the unpolar- 
ized differential cross section we need the sum the above 
terms and for the polarized one we need the difference 
similar to the three gluon case. Since the sum is over 
over all the polarization states we have to multiply by 
2 to include the CP conjugates. Then divide by 32 to 
average over the initial gluon colors and initial gluon and 
photon polarizations. Finally we have to divide by 167rs 2 . 



1. Color Singlet 

For 1 S , 3 P , 3 Pi and 3 P 2 we find 

|Af(+,+,+)| 2 = |A/(+,+,-)| 2 = 
\M (+, -, +)| 2 = \M (-, +, +)| 2 = . (49a) 
For 3 Si we find 

|A/(+,+,+)| 2 = (50a) 
128. 9 4 e 2 (i? [ 3 S[ 1] ]) 



M 2 s 2 (t + u) 2 
X (Q - Af 2 P) 2 

m2 128 5 4 e 2 <P[ 3 ^i 1) ]^ 



3ttM 
M 2 u 2 {s + t) 2 
X (Q - M 2 P) 2 

128g 4 e 2 (J?[ 3 5^ ) ]; 
3ttM 
M 2 t 2 (u + s) 2 
X (Q - M 2 P) 2 ' 



(50b) 



(50c) 



(50d) 



Our results agree with the upolarized and polarized re- 
sults in (A. 3) in KKMS if we make the replacement 

{R[ 3 S{ 0) }} = 16vr(0 [ 3 Si 0) ])/3. 
For 1 Pi we find 



|Af(+,+,+)| 2 = 



2 _ 1024 ff 4 e 2 (i?[ 1 P 1 (1) ] 



7tA/ 3 



Af 10 (5Q - Af 2 P) 



\M(+,+,-)\ 



(Q - APP) 3 

2 _ 1024 3 4 e 2 (P[ 1 P 1 (1) ]) Af 2 s 2 



(51a) 



tM 3 



(q - Nppy 



: [5A/ 4 Q - 4PQ - A/ 6 P - 2Q(i 2 + u 2 )] (51b) 



1024< 7 4 e 2 (P[ 1 P 1 (1) ]) M 2 



(q - Nppy 



x [5M 4 Q - 4PQ - Af 6 P - 2Q(s 2 + t 2 )] (51c) 



, |2 1024 5 4 e 2 (P[ 1 P 1 (1) ] 
|Af(-,+,+)| 2 = y 



APt 2 



ttA/ 3 (Q - APP) 3 

x [5A/ 4 Q - 4PQ - APP - 2Q(s 2 + u 2 )} . (51d) 

When we calculate the differential cross sections they 
agree with the unpolarized and polarized results in (A. 4) 
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in KKMS if we make the replacement (R^P^]) = 
8t:(0 pif ] }/9. 



2. Color Octet 



For 1 5'o we find 



2 _ 96g 4 e 2 (R[ 1 S^ ) }) usM 8 



iM(+,-,+)r 



7rM 


t(Q 


- M 2 P) 2 

(52a) 






4 

uss 


7rM 


i(Q 


- M 2 P) 2 

(52b) 






usu 4 


7rM 


i(Q 


- M 2 P) 2 

(52c) 








7rM 


t(Q 


- M 2 P) 2 ' 



(52d) 

When we calculate the unpolarized differential cross sec- 
tion, it does not agree with the first terms in (A. 5) in 
KKMS. However the polarized differential cross section 
agrees with the second terms in (A. 5) in KKMS if we 

make the replacement (R^S^}) = 2tt(0[ 1 S^ ) }). The 
reason for the different results is again caused by the fact 
that ghost contributions were not included in the KKMS 
analytic answers but are included in the KKMS fortran 
programs. Both the sum and the difference agree with 
(Af) and (A2) in YDHC. The sum agrees with (Al) in 
KLS. 

For 3 Si we find 

|M(+,+,+)| 2 = 



\M(+,+, 



|M(+,-,+)| 2 = 



iM(-,+,+)r 






(53a) 


320 5 4 e 2 (P [ 3 S{ 8) ] 


) s 2 M 2 (t + u) 2 


3ttM 


(Q-M 2 P) 2 




(53b) 


320g i e 2 {R [ 3 s{ 8) ] 


)u 2 M 2 {t + s) 2 


3ttM 


(Q-M 2 P) 2 




(53c) 


320 ff 4 e 2 (P [ 3 s{ a) ] 


)t 2 M 2 {s + u) 2 


3ttM 


(Q - M 2 P) 2 ' 




(53d) 



The sum and the difference both agree with (A. 6) 
in KKMS if wc make the replacement (R [ 3 Sf } ] ) = 

2 7 r(0[ 3 sf ) ])/3. They also agree with (A4) and (A5) 
in YDHC. 



2560 g 4 e 2 (i? pPf)]) 
irM 3 
M 10 (5Q - M 2 P) 
X (Q - M 2 P) 3 

2560 3 4 e 2 (P[ 1 P 1 (8) ]) M 2 s 2 



(54a) 



For 1 Pi we find 



\M(+,+,+)\ 2 = 



\M(+,+,-)\ 2 = 



|M(+,-,+)| 2 = 



|M(-,+,+)| 2 = 



The sum and the difference both agree with (A. 7) 
in KKMS if wc make the replacement (i?[ 1 P 1 (8) ]) = 

n(0M 6) ])/9- 
For 3 P we find 

2 l28g^ 2 (R[ 3 P^])9usM%M 2 -s) 2 

\M(+,+,+)\ ^3 t (Q-M 2 P) 4 

t 2 u 2 + 2Q(M 2 - s) + s 2 (M 2 - s) 2 





irM 3 (Q - 


M 2 P) 3 


X 


3M 4 Q - M 6 P + 2Qs 2 


(54b) 


2560 5 4 e 2 (^[ 1 ^i 8) ]> M 2 u 2 




ttM 3 (Q - 


M 2 P) 3 


X 


3M 4 Q- M 6 P + 2Qu 2 


(54c) 


2560 5 4 e 2 (P4 1 ^i (8) ]} M 2 t 2 




ttM 3 (Q - 


M 2 P) 3 


X 


3M 4 Q-M 6 P + 2Qi 2 


. (54d) 



+2sQ + 2s 3 {M 2 -s 



(55a) 



\M(+,+ 



2 128g 4 e 2 (P [ 3 P (8) ] ) uss 2 (M 2 - s) 2 
ttM 3 t(Q - M 2 P) 4 

At 2 u 2 {M 2 - sf + 4Qtu(M 2 - s) + Q 2 



-UsQ(M 2 - s) 2 - l8s 2 Q(M 2 - s) + 9s 6 
+9s 4 (M 2 - s) 2 - 6s 3 Q + 18s 5 (M 2 - s)l (55b) 

2 128g 4 e 2 (P[ 3 P (8) ]) u S u 2 (M 2 -u) 2 
1 l+ ' ' +JI '" 7rM 3 t(Q-M 2 P) 4 

At 2 s 2 {M 2 - u) 2 + 4Qts(M 2 -u) + Q 2 

-12uQ(M 2 - uf - 18u 2 Q(M 2 - u) + 9m 6 

+9m 4 (M 2 - u) 2 - Qu 3 Q + 18w 5 (M 2 - u)] (55c) 

| M( _ +j+) |2 = 128. 9 4 e 2 (P[ 3 P (8) ])^ 4 (M 2 -t) 2 



ttM 3 t(Q - M 2 P) 4 

4(M 2 - t) 4 + Asu(M 2 - t) 2 + s 2 u 2 
+28t(M 2 - t) 3 + UQ(M 2 -t) + 25t 4 
+69t 2 (M 2 - i) 2 + lOiQ + m 3 {M 2 - i) 



(55d) 



and only the difference agrees with (A. 8) in KKMS if 
we make the replacement (R [ 3 P (8) ] ) = 2tt(0 [ 3 P (8) ] ). 
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The sum and the difference agree with (A6) and (A7) in 
YDHC once a typo is corrected; the last term in these 
equations should have been (t + u)~ 2 instead of (t + s)~ 2 . 
The sum agrees with (A2) in KLS. 

For 3 Pi we find 



|M(+,+,+)| 2 = (56a) 

2 _ 1152/e 2 (P[ 3 P^])s 2 (M 2 -s) 2 
1 >l ~~ ttM 3 (Q-M 2 P) 4 

x s 5 (t + u) 2 + s{9t 2 u 4 + 20tV + 13iV + 2t 5 u) 

+s 2 (tu 4 + 26iV + 38t 3 u 2 + I3t 4 u + t 5 + u 5 ) 
+s 3 (6tu 3 + Mt 2 u 2 + 22t 3 u + 3t 4 - u 4 ) 

+s 4 (9tu 2 + I3t 2 u + 3t 3 - u 3 ) + t 2 u 2 (t + u) 3 ] (56b) 

2 _ ll52g 4 e 2 (R[ 3 Pl S) ])u 2 (M 2 -u) 2 

I l+ ' ttM 3 (Q-M 2 P) 4 

x u 5 (s + t) 2 + u(2st 5 + 13s 2 i 4 + 20s 3 t 3 + 9s 4 t 2 ) 

+m 2 (13s< 4 + 38s 2 i 3 + 26s 3 t 2 + s 4 t + s 5 + t 5 ) 
+u 3 {22st 3 + 34s 2 t 2 + 6s 3 t - s 4 + M 4 ) 

+w 4 (13st 2 + 9s 2 t - s 3 + 3i 3 ) + s 2 t 2 (t + s) 3 (56c) 

2 1152g 4 e 2 (R[ 3 P^})t 2 (M 2 -t) 2 

II ' + ' +JI ~~ ttM 3 {Q-M 2 Py 

x ts 2 u 2 (5(u + s) 2 - 2su) - 4s 3 u 3 (u + s) 
+t 2 ((u + s) 5 + 8s 2 u 2 (u + s)) + s 2 u 2 (u + s) 3 



+t 3 (3(u + s) 4 - 2su(s 2 + u 2 )) + t 5 (s 2 + u 2 ) 
+t 4 (3(u + s) 3 - Asu{u + s)) 



(56d) 



and only the difference agrees with (A. 9) in KKMS if 

we make the replacement (R[ 3 P^ S) }) = n(0 ppf 5 ] )/4. 
The sum and the difference agree with (A8) and (A9) 
respectively in YDHC. The sum also agrees with (A3) in 
KLS once the factor (s 2 — u 2 ) 2 is replaced by (s 2 — u 2 ) 4 . 

For 3 P 2 we find that only the difference of the results 
below agrees with (A. 10) in KKMS if we make the re- 
placement (i?[ 3 P 2 (8) ]) = 16tt(0[ 3 P 2 (8) ])/15. The differ- 
ence also agrees with (A10) and (All) in YDHC. The 
sum from our results does not agree agree with the an- 
swer in (A4) in KLS as well as with (A10) in YDHC. 
The diferent results in (A.5), (A.8), (A.9) and (A.10) of 
the KKMS paper are again due to the inadvertent omis- 
sion of gluon contributions to their answers. They say 
they have included these contributions in their fortran 



programs. 

|M(+,+,+)| 2 = (57a) 

, _. |2 48. 9 4 e WP, (8) ]) s 2 u(M 2 -s) 2 
1 >l ' 7rM 3 Q(Q-M 2 Pf 

x (12*V + 48t 3 u 6 + 72t 4 u 5 + 48t 5 u 4 + 12i 6 u 3 ) 

+s{2Atu 7 + 96t 2 u 6 + I7lt 3 u 5 + mt 4 u 4 + 105tV 
+27* V) + s 2 (72tu 6 + 140* V + 187* 3 w 4 + 200* V 
+111*V + 18* 6 m + 12m 7 ) + s 3 (51*m 5 + 59*V 
+ 134* 3 w 3 + 162* 4 u 2 + 63* 5 u + 3* 6 + 24m 6 ) 
+s 4 (-3*u 4 + 26* V + 102t 3 u 2 + m 4 u + 9t 5 
+ I2u 5 ) + s 5 (-3*u 3 + 27*V + 39t 3 u + 9t 4 ) 

(57b) 



+s 6 (3*w 2 + 6t 2 u + 3* 3 ) 
|M(+,-,+)| 2 = 



48g 4 e 2 {R [ 3 P 2 (8) ] ) su 2 (M 2 - u) 2 



ttM 3 Q(Q - M 2 P) 4 

(12s 3 * 6 + 48s 4 * 5 + 72s 5 * 4 + 48s 6 * 3 + 12s 7 * 2 ) 

+w(27s 2 * 6 + 105s 3 * 5 + 177s 4 * 4 + 171s 5 * 3 + 96s 6 * 2 
+24s 7 *) + u 2 (18s* 6 + Ills 2 * 5 + 200s 3 * 4 + 187s 4 * 3 
+ 140s 5 * 2 + 72s 6 * + 12s 7 ) + u 3 (tt 6 + 63si 5 + 162s 2 i 4 
+ 134s 3 t 3 + 59s 4 t 2 + 51s 5 t + 24s 6 ) + u 4 (9t 5 + 78st 4 
+ 102s 2 t 3 + 26s 3 t 2 - 3s 4 t + 12s 5 ) + u 5 (9t 4 + 39st 3 

+27 s 2 t 2 - 3s 3 1) + w 6 (3t 3 + 6si 2 + 3s 2 t)} (57c) 



\M(-,+,+)\ 2 



48 5 4 e 2 (P[ 3 P 2 (8) ]) su{M 2 -tf 



ttM 3 Q(Q - M 2 P) 4 

x 24Qsu(u 5 + 5su 4 + 10s V + 10sV + 5s 4 m + s 4 ) 

+ l2Qt(u 6 + 10su 5 + 29sV + 40s 3 w 3 + 29sV 
+ 10s 5 u + s 6 ) + 3Qi 2 (16u 5 + 89su 4 + 183s 2 m 3 
+ 183sV + 89s 4 m + 16s 5 ) + Qt 3 (92u 4 + 367su 3 
+552s 2 u 2 + 367s 3 m + 92s 4 ) + t 5 (3u 5 + 119su 4 
+334s 2 u 3 + 334s 3 w 2 + 119s 4 u + 3s 5 ) + t 6 (9u 4 
+ 102su 3 + 182s 2 u 2 + 102s 3 w + 9s 4 ) + t 7 (9u 3 
+47su 2 + 47s 2 w + 9s 3 ) + t 8 (3u 2 + 8su + 3s 2 ) 



+12sV(it 4 + 4stt 3 + 6sV + 4s 3 u + s 



2„.2 



(57d) 



IV. CONCLUSIONS 

We have calculated the gluon-gluon and photon-gluon 
amplitudes for the production of color singlet and color 
octet charmonium production. These amplitudes are re- 
quired for the QCD analysis of charmonium production 
in polarized and unpolarized hadron-hadron and photon- 
hadron collisions. Our calculations clarify several incon- 
sistencies in previously published results. 
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